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Special solutions of the Hartree-Fock (HF) problem for Coulomb interacting electrons, being 
described by a simple model of the Cu-O planes in La2Cu04, are presented. One of the mean field 
states obtained, is able to predict some of the basic properties of this material, such as its insulator 
character and the antiferromagnetic order. The natural appearance of pseudogaps in some states 
| of this compound is also indicated by another of the HF states obtained. These surprising results 

follow after eliminating spin and crystal symmetry restrictions which are usually imposed on the 
1 single particle HF orbitals, by means of employing a rotational invariant formulation of the HF 

scheme which was originally introduced by Dirac. Therefore, it is exemplified how, up to now being 
considered strong correlation effects, can be described by improving the HF solution of the physical 
■ systems. In other words, defining the correlation effects as such ones shown by the physical system 

£N) ' and which are not predicted by the best HF (lowest energy) solution, allows to explain currently 

assumed as strong correlation properties, as simple mean field ones. The discussion also helps to 
clarify the role of the antiferromagnetism and pseudogaps in the physics of the HTSC materials and 
D ' indicates a promising way to start conciliating the Mott and Slater pictures for the description of 

^ l the transition metal oxides and other strongly correlated electron systems. 

> . 

^ 1 PACS numbers: 71.10.Fd,71.15.Mb,71.27.+a,71.30.+h,74.20.-z,74.25.Ha, 74.25. Jb,74.72.-h 
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The Hubbard models in the theory of strongly correlated electron systems are notably successful [E S S & 
Hi S H @, They remarkable reproduce the properties of Mott insulators, such as metal-transitions oxides and 
copper-oxygen layered HTc compounds Q. However, the search for approaches having more basic foundations do 
not cease, hoping that they could open the way for obtaining more exact and specific results ll|. The so called 
first principle or Slater procedures, are electronic structure calculations that begin by considering the interactions 
among electrons or atoms in vacuum. The study of the band structure they predict, in principle should offer a way 
£%J . for the precise determination of the physical properties of each material [ll , T^i ll[ • However, the above mentioned 
^ ' potentialities had been failing in describing many systems showing strong correlation effects [Tj. 

The motivation of the present letter emerged from a suspicion that the Hartree Fock (HF) method, could had been 
underestimated in its possibilities for helping in clarifying the above mentioned difficulties @, [3] • As a net result of 
this work, we came to believe that a large deal of correlation effects, can be described in the framework of the HF 
scheme, after removing certain symmetry restrictions which obstacle the finding of the best HF solutions. By example, 
it has been early shown by Slater in Ref. |2j, that sometimes the HF potential breaks the symmetry of the original 
qq ' crystalline lattice, creating magnetic structures and gaps. This symmetry breaking effect has been also more recently 
underlined and deepened in Ref. [l6| . However, the removal of the lattice symmetry restrictions alone had not been 
able to describe the insulator properties of a large class of materials [2, 0, E| ■ One of the central results of the present 
investigation, is the identification of another kind of symmetry restriction that seemingly had been overlooked for a 
long time. It corresponds to the usual assumption about the necessary a (S z = 1/2) or /3 (S z = —1/2) orientations 



of the spin of each solution for the HF orbitals [12J, |17|. For supporting this main statement, the work considered 
the HF problem as applied to a simple one band model of the superconducting material La-aCuO,! [15l|. looking from 
the start for single particle orbitals being non separable in their spacial and spinor dependence. That is, they will 
have the structure <fr(x, s) ^ 4>{x)^{s), i.e. the orbitals have not an absolute common quantization direction for their 
electron spin. Note, that to proceed in this way means not other thing that to apply the Dirac's formulation of the 
HF procedure [l3 |. The results, as it will be seen, are able to dscribe the basic properties of the La2Cu04 (lij . 
We would like to stress that the obtained results do not look so radical if we consider the following circumstance: 
the correlation effects are associated with the difference shown by the system's properties with the ones predicted 
by the HF procedure applied to it. Therefore, what we have argued here, is simply that the there exist unexplored 
improvements of the mean field schemes which allow to predict properties which are currently considered as strong 
correlation properties. These modifications of the HF method are related with usual constraint which are imposed 
on the crystalline and spin properties of the single particle orbitals. The possibility of their employment for solving 
the longstanding debate between the Slater and Mott pictures in the theory of transition metal oxides and HTc 



2 



superconducting materials will be considered elsewhere. 

In this letter, we will firstly describe the physical basis of the model employed and after that, the results obtained 
for the various HF solutions will be presented. For the details of the concrete solutions we refer to the extended 
version of this work in Ref . [3| ■ 

Let H(xi, ...xn) = X)i^o(^i) + \ Tlj^i V{ x ii x j)> be the N-electrons system hamiltonian, including kinetic plus 
interaction with the environment hamiltonian H.q, besides Coulomb interaction among pairs of electrons V. The HF 
equations of motion were solved by allowing for nonseparable solutions for the orbitals 4>r){x, s), where rj — ki, Un 
are the labels of the HF basis they form. The system of HF equations is rotational invariant when formulated in the 
general HF procedure firstly introduced by Dirac in Ref. (13 ] . Its explicit expression for the here considered problem 
can be found in Ref. [II [01. 

Let us first present the effective band model used to describe the dynamic of the less bounded electrons in La2CuC>4. 



It is known that at low temperature La2CuC>4 is an antiferromagnetic-insulator 15j. However, in evident contradiction 



with the experiments the Linear Augmented Plane Waves (LAPW) band calculations predicts for it metal and para- 
magnetic zero temperature properties 12|. Nevertheless, such band studies results at least show that the conduction 
electrons are strongly coupled to the Bravais lattice centers of the copper oxygen planes. Clearly this tight-binding be- 
havior is determined by the interaction of the electrons with its surrounding effective environment. The less bounded 
electron in the La2CuC>4 molecule is the Cu 2+ 's not paired one. At difference from O 2- ions, Cu ones do not have 
their last shell (3d) closed. Those copper 3d electrons partially fill the last band of La2Cu04 solid and in what follows 
they shall be referred as: the electron gas. It seems appropriate to consider those electrons as strongly linked to 
Cu02 cells and with special preference for the Cu centers Thus, our Bravais lattice is going to be the squared net 
coincident with the array of copper sites (see the figure[T]). Further, the presence of electrons pertaining to the various 
fully filled bands in the material plus the nuclear charges, will play a double role in the model. Firstly, they will act 
as an effective polarizable environment screening the field of the charges in the electron gas. It will be reflected by a 
dielectric constant e reducing the Coulomb interaction. Secondly, the mean field created by the environment will act 
as a periodic potential W 7 , tight-binding the electrons to the Cu centers. The interaction F b among the electron gas 
and a "jellium" neutralizing its charges is also considered. It was be modeled as a gaussian distribution of positive 
charges Pb(y) = ^p- exp(— ^-), surrounding each lattice point and having a characteristic radius b . In resume, the 
free hamiltonian of the model takes the form 
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H (x)=^i^ + ^ 7 (x) + F b (x), (1) 

i=l 

VF 7 (x) = W 7 (x + R), 

,2 



F bW = V / d y— — , b < p, 

4iree ^J x-y 



R 

where pf is the i-th electron's squared momentum operator; m is the electron mass; eo is the vacuum permitivity and 
lattice vectors R = n Xl pe Xl -\- n x . 2 pe x . 2 with n Xl and n X2 being integers, move on Bravais lattice. It will be referred 
in what follows as the the absolute lattice. The versors e Xl and e X2 are resting on the direction defined by the lattice's 
nearest neighbors (see figure [I] a)). The distance between Cu nearest neighbours is p ~ 3.8 A [lj|. The interaction 
among pairs of electrons is taken in the form V(x, y) = 4 ^ eeo | x -y| ' wm ch includes the dielectric constant associated 
to the presence of the effective environment. 

We are seeking here for HF solutions with orbitals having a non separable spin and orbit structures. Thus, it 
was considered that the spin can show a different projection for the different Wannier orbitals to be superposed for 
defining those orbitals. The spin for each of them will be either or a or (3 type, according they are linked either to 
one or the other of the two sublattices shown in figure [T] a). The points of these sublattices were defined as follows: 
R( r ) — \/2nip qi + \f2ri2P q2 + q*- r \ r = 1, 2; with m and U2 S Z , and in which the vector q( r ) = if r = 1 
and q( r ) = p e Xl when r = 2; and where qi and q2 form the basis versors on each one of them. The searched 
solutions were constructed as eigenfunctions of the operators T R ( r ) belonging to the translation group transforming 
each sublattice on itself: T R ( r )(/)k,; = exp(i k • R( r ))0k,i- Therefore we imposed periodic boundary conditions on the 
4>k,i in the absolute lattice's boundaries xi= - p L and p L, L , X2= -p L and p L (see figure [T] a) ) . This condition 
determines the allowed set of quasimomenta k =fp + n 2&x 2 )i with Ux, n 2 G Z and — ^ < ni ± < ^. Note 

that we are demanding less than the full crystal symmetry on the single particle states which we are looking for. Let 
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FIG. 1: The figures show: a) The point lattice associated to the Cu-O planes. For removing the symmetry restrictions, it will 
be helpful to separate the lattice in the two represented sublattices; and b) shows the corresponding base of the Cu-O planes. 



the single particle states represented in the explicitly non separable form 

k , i (x, s ) = E B ^!^ ) ( x ' s )' 



where I is the additional label needed for indexing the stationary states. The tight binding Bloch basis tp^ 
appearing is defined as 



x, s) = u°* (s) J2 cxp(* k • R«) few (x), 

RW 



cr,u * — a-, u 



Vrm ( x ) 
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where N is the number of electrons in the electron gas, a z is the spin z projection operator, where z is the orthogonal 
direction to the copper oxygen Cu02 planes; o~ z = -1, 1, are the eigenvalue of the previously mentioned operator and 
r = 1, 2, is the label which indicates each one of the sublattices. As we are working at half filling condition, then N is 
equal to the number of cells in the periodicity region N c . Note that int his simple model, the electron were considered 
as moving in two dimensions and that gaussian T.B. Wannier orbitals had been assumed [HI]. 
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FIG. 2: The figure 1 a) shows the Brillouin zone associated the absolute lattice. The grey region signals the Brillouin zone 
(B.Z.)of the sublattices. The unit of quasimommentum is ^. Figure b) shows the doubly degenerated bands associated to 
a paramagnetic and metallic state. The zero energy level in all the band diagrams is the Fermi energy of the isolator and 
antiferromagnetic solution. 

The first solution considered, to be called the PM one, was one having their orbitals being eigenfunction of the 
maximal group of translations leaving invariant the absolute lattice. That is, the maximum possible crystal symmetry 
is demanded. Also, the spin structure of the orbitals were assumed to be of a or (3 types. For this case, figure[2]shows 
the paramagnetic, metallic and doubly degenerate band obtained from an iterative process of solving the HF equation 
in Ref. [l^|. For the case of N = 20 x 20 electrons, the momenta of occupied states are shown in figure [5] a). They lay 
inside the shadowed square of side V2tt/p: the Brillouin zone (B.Z) of the sublattices. The condition of reproducing 
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the bandwidth of 3.8 eV for the conduction band reported in Ref. [l2[, allowed to fix the dimensionless values of 
the free parameters of our model to be: e=10, a=0.25, &=0.05 and 7=-0.03. See Ref. [HI for those dimensionless 



definitions. The here obtained band topologically coincides with the conduction band result of Ref. [12 1. 
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FIG. 3: Energy bands obtained for: a) A sample of 20x20 cells, Eg ap = 1.32 eV. b) A sample of 30x30 cells, E gap = 1-32 eV. 
The zero energy level is chosen in the Fermi level of the 20x20 system and the plotting region is the B.Z. of the sublattices: 
the shadowed square of side \/2n/p in Fig. 2 a). 



After being determined the parameters of the model, we searched for general HF solutions showing a non separable 
structure, and moreover, also retaining a reduced translational symmetry in the defined sublattices. The iterative 
process employed for fin ding this solution, started from a particular state having an antiferromagnetic character from 
the beginning (see Ref. [18|). This initial form helped the convergence toward the results being depicted in figure [3] 
It shows two sets of bands following for a half filling band, associated to two periodicity lattices having 20x20 and 
30x30 cells. The difference between the energies in them is of the order of 10 dimensionless units of energy JL ^ = 
8.3 eV. The measure of the gap shown indicate that the bands correspond to insulating states. The close similarity 
between both results indicates that the thermodynamical limit has being achieved. The HF energy of this solution 
was the lowest among of all the ones found. It also follows that in coincidence with the experimental evidence, these 
states show a local magnetic moment resting on the direction of the sublattice X12 as shown in figure 2J In the figure 
[4] a) the only non vanishing component of the magnetization m in this solution is plotted. An interesting result 
is that it has been experimentally observed that La2CuC>4 has a magnetic moment of 0.68 /is per Cu site on the 
CuO plane [loj ]. and the value obtained from evaluating the magnetization of this HF state turns out to be 0.67 /is. 
Therefore, the considered here discussion satisfactorily predicts the whole antiferromagnetic structure of La2CuC>4. 
The antiferromagnetic and isolator solution will be referred below as the IAF state. The single particle states of 
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FIG. 4: The magnetization vector m of the more stable HF state rests in the direction 1-2. a) The figure shows the projection 
fni2 of the dimensionless magnetization, in the 1-2 direction. The magnetization unit is ^f-. b) The picture shows a scheme 
of the mean magnetic moment per site in the lattice. The figure insertion below depicts the angle between the magnetization 
components of each orbital on each of the sublattices. Note that the orbitals are more AF like as they become close to the 
Fermi level. The plotting area is again the B.Z. of the sublattices shown in Fig. 2 a). 

this solution carry a more intensive antiferromagnetism as more closer they are in energy from the Fermi surface. 
This property offers a clear explanation of the gradual loss observed in the antiferromagnetic order under the doping 
with holes [7| . The dependence of the angle <j> between the magnetic moments per cell on each of the two sublattice 
shown by a given Bloch orbital, is plotted in the insertion at the bottom of figure 2] These components are defined 
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as the integrals of the magnetic moment over all the unit cells of the absolute lattice centered in the sublattice points 
[ill]. Then, this HF solution indicates that after the orbitals are allowed to spatial dependent spin orientations, the 
electrons prefer to reorient their spin when traveling between contiguous lattice cells. Note that the states laying just 
on the Fermi surface are perfectly antifcrromagnetic ones, and that the more away from the boundaries the orbitals 
are, the less antiferromagnetic they become. Further, in figure a) the band spectra corresponding to another HF 
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FIG. 5: a) The band structure associated to the paramagnetic ground state also showing a pseudogap. b) The figure shows 
in the same plot the two occupied bands linked with the pseudogap and the IAF states. The difference in their energies is 
concentrated near the Fermi level. The zero energy reference coincides with the Fermi level of the IAF state in both pictures. 
The domains also in both cases is the B.Z. of the sublattices shown in Fig. 2 a). 

paramagnetic ground state (to be named as the PPG one), obtained from the HF equations in Ref. [18[ is shown. This 
solution was obtained by only requiring the orbitals to show the full translational symmetry of the CuO planes but 
not an a or (3 type of their spin structure. Note that simply allowing this last freedom determined the existence of a 
pseudogap in the state. It reaches a maximum value of 0.1 eV « 10 3 K (equivalent to 0.012 dimensionless unit of energy 
^2 = 8.3 eV). It is an interesting outcome that the HF energy of this ground state is exactly coincident with the one 
corresponding to the paramagnetic and metallic state. Moreover, the occupied single particle states in both solutions 
arc identical and in consequence the momentum dependence of the filled energy bands also coincide. Henceforth, the 
difference between the two solutions only refers to the non occupied states. It thus follow that the HF energies per 
particle of the paramagnetic- metallic (PM), and the paramagnetic with pseudogap (PPG) HF states coincide. They 
show a value +0.0761 eV higher than the energy per particle corresponding to the insulator-antiferromagnetic (IAF) 
ground state. The spin structure of the excited states can be seen as a dynamical "correlation" effect determined by 
the freedom allowed to the HF orbitals. It is worth noticing that the energy difference PM (PPG)-IAF and the Neel 
temperature of this kind of materials are both of the order of 10 2 K. Then, the results suggest the possibility of having 
further success in applying the approach started in this work to the description of some regions of the phase diagram 
of the La2 CUO4. In Figure Ob) the PPG (PM) and IAF occupied bands are depicted in a common frame. The main 
difference in their energies corresponds to the single particle states being closer to Fermi surface. As noted before the 
same behavior has the antiferromagnetic character of the single particle states of the IAF solution. Therefore, it is 
suggested that both solutions could evolve toward a common ground state lacking an absolute magnetic order under 
under doping with holes. This indication is in correspondence with the pattern shown by the phase diagrams of these 
materials [la]. 
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